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We study exotic fermions with spectrum E ot p . Such spectrum emerges in the vicinity of the Fermi 
point with multiple topological charge N, if special symmetry is obeyed. When this symmetry is violated, the 
multiple Fermi point typically splits into N elementary Fermi points - Dirac points with N = 1 and spectrum 
E 2 oc p 2 . 



PACS: 

1. INTRODUCTION 

There is a fundamental interplay of symmetry and 
topology in physics, both in condensed matter and rel- 
ativistic quantum fields. Traditionally the first role was 
played by symmetry (symmetry classification of crys- 
tals, liquid crystals, magnets, superconductors, super- 
fluids, etc.). The phenomenon of spontaneously broken 
symmetry remains one of the major tools in physics. 
The last decades demonstrated the opposite tendency 
in which topology is primary. Topology in momen- 
tum space is becoming the main characteristics of quan- 
tum vacua - ground states of the system at T = 
HI ® 1 a E3 G3 E M M HU El HI D3- For ex- 
ample in topological matter, nodes in the spectrum of 
fermionic quasiparticles are protected by topology, and 
this property is insensitive to the details of the micro- 
scopic many-body Hamiltonian (see e.g. review [5]). 
The momentum-space topological invariants determine 
universality classes of the topological matter and the 
type of the effective theory which emerges at low energy 
and low temperature. In many cases they also give rise 
to emergent symmetry. Examples are provided by the 
nodes of the energy spectrum characterized by the ele- 
mentary topological charges, N — +1 or N = — 1. Close 
to such nodes the effective Lorentz invariance emerges: 
the Fermi point with N = +1 or N — — 1 represents 
the Dirac point and the spectrum forms the relativistic 
Dirac cone. This is the consequence of the Atiyah-Bott- 
Shapiro construction [T3] . 

However, in many systems (including condensed 
matter and relativistic quantum vacua), the Fermi 
points with N = +1 or N = — 1 may merge together 
forming the points either with N = or with multiple 



N (i.e. | TV | > 1 [H]). In this case topology and symme- 
try become equally important, because it is the symme- 
try which may stabilize the degenerate node. Example 
is provided by the Standard Model of particle physics, 
where 16 fermions of one generation have degenerate 
Dirac point at p = with the trivial total topologi- 
cal charge iV = 8 — 8 = 0. In the symmetric phase 
of Standard Model the nodes in the spectrum survive 
due to a discrete symmetry between the fermions and 
they disappear in the non-symmetric phase forming the 
fully gapped vacuum [8]. In case of degenerate Fermi 
point with \N\ > 1, situation is more diverse. Depend- 
ing on symmetry, interaction between fermionic flavors 
may lead to splitting of the multiple Fermi point to ele- 
mentary Dirac points |17) ; or gives rise to the essentially 
non-relativistic energy spectrum E±(p — > 0) — > ±p N , 
which corresponds to different scaling for space and time 
in the infrared: r — > Ar, t — > X N t. The particular case 
of anisotropic scaling with N = 3 was suggested by Ho- 
rava for quantum gravity at short distances [TBI HI HO] > 
and anisotropic scaling in the infrared in Ref. |21j . 

The non-linear spectrum arising near the Fermi 
point with N — 2 has been discussed for different sys- 
tems including graphene, double cuprate layer in high- 
T c superconductors, surface states of topological insu- 
lators and neutrino physics [521 18] [23 1 1^ 1 I^ H 
[551[55]. The spectrum of (quasi)particles in the vicin- 
ity of the doubly degenerate node (say, with topological 
charge N = N% = ±2 for Fermi points of co-dimension 
3) depends on symmetry. One may obtain: two Weyl 
fermions, if there is some special symmetry; exotic mass- 
less fermions with quadratic dispersion at low energy, 
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or semi-Dirac fermions with linear dispersion in one di- 
rection and quadratic dispersion in the other, 



£ ± ( P )«±Wc 2 p2 + 



Pi 

2m 



(2) 



Here we discuss the class of effective Hamiltonians 
which has a Fermi point with higher degeneracy, de- 
scribed by the symmetry protected topological invari- 
ant N — 3 and N = 4, and quantum phase transitions 
which may occur in these systems where topology of the 
spectrum changes. 



2. CUBIC SPECTRUM 

Let us consider first the case with N = 3. Examples 
are three families of right-handed Weyl 2-component 
fermions in particle physics; three cuprate layers in high- 
T c superconductors; three graphene layers, etc. If the 
Fermi point is topologically protected, i.e. there is a 
conserved topological invariant N, the node in the spec- 
trum cannot disappear even in the presence of inter- 
action, but it can split into N nodes with elementary 
charge N = 1. The splitting can be prevented if there is 
a symmetry in play, such as rotational symmetry. Here 
we provide an example of such symmetry. 

For simplicity we study the 2+1 systems. In general 
the nodal points in 2+1 dimensions (Fermi points of co- 
dimension 2) obey Z 2 topology that is why we 
also need an additional symmetry K which extends the 
group Z 2 to the full group Z and thus makes the mul- 
tiple Fermi point possible. The corresponding invariant 
protected by symmetry K [231 130] is: 



N 



K 



— tr idl KG(uj = 0, p)d l G- 1 (ui = 0,p) , 
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(3) 

where G is Green's function matrix, and G is contour 
around the Fermi point in 2D momentum space (p x ,p y ), 
or around the Fermi surface if the Fermi point expands 
to the Fermi surface. The matrix K commutes or anti- 
commutes with matrix G{uj = 0,p). The single particle 
Green's function at zero energy represents the effective 
single-particle Hamiltonian H(p) = G _1 (cj = 0,p). 

We consider 3 species (families or flavors) of 
fermions, each of them being described by the invariant 
N = N K = +1 and an effective relativistic Hamiltonian 
emerging in the vicinity of the Fermi point 



H (P) = CT ■ P = (T X P X + (JyPy 



The matrix K = cr z anticommutes with the Hamilto- 
nian. This supports the topologically protected node in 
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Fig. 1 Spectrum of the Hamiltonian (J7]) showing cu- 
bic dispersion for the lowest two eigenvalues around 
the point p = 0. Different colors correspond to dif- 
ferent eigenvalues, spectrum is shown as function of p x 
at p y = 0. The spectra have been calculated with equal 
coupling strengths g 12 = gia = #23 = g- 

spectrum, which is robust to interactions. The position 
of the node here is chosen at p = 0: 



E 2 = p 2 . 



(5) 



The total topological charge of three nodes at p = of 
three fermionic species is N K = +3. Let us now intro- 
duce mixing of the fermions, which obeys some symme- 
tries which may follow from the underlying microscopic 
theory. We consider two cases 



m(p) 
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g 2 \o~ cr ■ p g 2 3<J + 
33ic + g 32 a~ cr ■ p 



(6) 



and 



n 2 (p) 
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where — \(o x ± i(j y ) are ladder operators. Both 
Hamiltonians anti-commute with K = a z and thus mix- 
ing preserves the topological charge N K in Hamil- 
tonian ([6]) is symmetric with respect to the Z3 group 
of permutations if all the couplings g mn are equal. For 
general couplings, Hamiltonian ^ is symmetric under 
the Z3 group of rotations by 27r/3 combined with gauge 
transformations from the SU(3) family group: 



^l(OorbP) = C/ g augc 

U sp i n rii{p)U spin U, 



gauge • 



(8) 



^ Here O or b and J7 sp i n are orbital rotations by 27r/3 in 



momentum and spin space respectively, and C/ gaug0 — 
diag (l, e 27 ™/ 3 , e " 2 ™/3) is the element of the S"/7(3) fam- 
ily group. 
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Fig. 2 Contour plot showing the equi-energy curves in 
the vicinity of the degenerate Fermi point with Nk = 3 
and cubic dispersion. Here gi2 = <?23 = g = (713/2. 

The Hamiltonian ([7]) at p x = p y = is independent 
of the spin rotations up to a global phase of the coupling 
constants. Under spin rotation by angle 9 all elements in 
the upper triangular matrix are multiplied by e 10 , while 
all elements in the lower triangular matrix are multi- 
plied by e~ te . This symmetry of triangular matrices 
generates specific property of the spectrum. Mixing be- 
tween fermions does not split the multiple Fermi point 
at p = 0, as a result the gapless branch of spectrum in 
Fig. Q]has the cubic form at low energy, E — >• 0, which 
corresponds to the topological charge N K — +3: 

£ 2 «7 3 V , 73 — I \ 1 ■ ( 9 ) 

1^12 1 1^23 I 

In the low-energy limit the spectrum in the vicinity of 
the multiple Fermi point (J9j) is symmetric under rota- 
tions. But in general the spectrum is not symmetric as 
demonstrated in Figs. [T]and[2J There is only the sym- 
metry with respect to reflection, (p x ,p y ) — > (p x ,—p y ). 
The rotational symmetry of spectrum Q is an emer- 
gent phenomenon. For trilaycr graphene this spectrum 
has been discussed in Ref. [31] . 

The cubic spectrum Q becomes singular when ei- 
ther one of the coupling constants 312 or (723 nullify. At 
this point of topological quantum phase transition, the 
gapless spectrum has two branches: one with quadratic 
dispersion E± = ±p 2 /2m and another one with linear 
dispersion E± = ±p. Note that in the limit p — >■ 
the spectrum ([9]) does not depend on (731, i.e. the spec- 
trum is not symmetric under permutations. This is the 



property of the Hamiltonian "^(p): there is no per- 
mutation symmetry even if the couplings are the same, 

.912 = 523 = 331- 



E/g 




Fig. 3 Spectrum of the permutation symmetric Hamil- 
tonian (|6]l shows three Dirac points (Fermi points with 
N = 1 and linear dispersion), located at {p x ,p y }/g = 
{-1,0}, {1/2, \/3/2},{l/2, -VS/2}, first of them is 
shown in the figure. Different colors correspond to dif- 
ferent eigenvalues, spectrum is shown as function of p x 
at p y = 0. The spectra have been calculated with equal 
coupling strengths 312 = 513 = #23 = g in (JTUJ . 
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Fig. 4 Splitting of Fermi point with iV = 3 into three 
Fermi points with N = 1 each in the spectrum of the 
Hamiltonian ©. We plotted the case of equal coupling 
strengths 312 = gis = #23 = g, corresponding to full 
permutation symmetry. 

At <73i = one also has the quantum phase tran- 
sition, since at 531 = the Hamiltonians (0 and ([6]) 
coincide: %2(p,33i = 0) = %i(p,<?3i = 0). Then 
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when (731 grows from zero in the multiple Fermi 
point splits into 3 elementary Dirac points, i.e. Fermi 
points with N K = +1. For equal coupling strengths 
3i2 = 913 = .923 = 3 one has 
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where p± = p x ± ip y . In this case the full permutation 
symmetry takes place and three Fermi points with N = 
1 each form the configuration obeying the three-fold ro- 
tational symmetry, see Figs. [3] and [4] For general cou- 
plings, zeroes are located at (p x + ip y ) 3 = -312323331, 
also obeying the Z3 symmetry in {5} . As one of the cou- 
plings vanishes, the Fermi points merge at p = and 
produce either a linear, quadratic or cubic dispersion 
relation, depending on which of the couplings vanish. 



FERMI POINTS WITH HIGHER 
DEGENERACY 
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Fig. 5 Spectrum of the Hamiltonian showing quar- 
tic dispersion for the lowest two eigenvalues around the 
point p = 0. Different colors correspond to different 
eigenvalues, spectrum is shown as function of p x at p y — 
0. The spectra have been calculated with equal coupling 
strengths g 12 = 313 = 323 = ffi4 = 324 = 334 = g- 



In case of four fermionic species, the mixing which 
does not produce splitting of the Fermi point is obtained 



by the same principle as in ([7]): all matrix elements 
above the main diagonal contain only tr + (or cr~): 
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Then again under spin rotation by angle 9 all elements in 
the upper triangular matrix are multiplied by e* e , while 
all elements in the lower triangular matrix are multi- 
plied by e~ 10 . As a result the multiple Fermi point with 
N = 4 is preserved giving rise to the quartic spectrum 
in vicinity of the Fermi point: 

1 
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74 



(12) 
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For the tetralayer graphene this spectrum was suggested 
in Ref. |32j . Again the rotational symmetry emerges 
new the Fermi point, but spectrum is not symmetric 
under permutations. In fact, there is no such 4x4 ma- 
trix that would consist of couplings described by ladder 
operators and which would be symmetric under permu- 
tations. 




Fig. 6 Splitting of Fermi point with N = 4 into four 
Fermi points with N = 1 each in the spectrum of the 
Hamiltonian where in Eq. (j 1 1 p the element proportional 
to the coupling 324 has been changed to o~ instead of 
tj + . All the (finite) couplings are assumed to have the 
same magnitude 3. 

In general the Fermi point with arbitrary N may 
jive rise to the spectrum 

E 2 



(13) 
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For graphene this spectrum is proposed in Ref. |33j . 
The discussed symmetry of matrix elements g mn ex- 
tended to 2N x 2N matrix gives (ITS")) with the prefactor 



MV I II | I | ■ 

|#12||ff23| ■ • ■ \9N-1,N\ 

Violation of this symmetry may lead to splitting of the 
multiple Fermi point into N elementary Fermi points 
- Dirac points with N = 1 and 'relativistic' spectrum 
E 2 oc p 2 . For N = 4 the splitting is shown in Fig. [51 

The effective Hamiltonian describing fermions in the 
vicinity of multiple Fermi point is H = <J~p+ + <r + p^ , 
see [21] • An example of the effective Hamiltonian 
describing the multiple Fermi point with topological 
charge N = N3 in 3+1 systems is [16] 

H = o z p z + o-pl + a+p N . (15) 

This Hamiltonian has the spectrum E 2 = p 2 + p± , 
which has linear dispersion in one direction and non- 
linear dispersion in the others. 

4. CONCLUSION 

Using the momentum-space topology we found the 
classes of Hamiltonians which exhibit the nonlinear 
Dirac- like spectrum E± = ±j]yp N . Such spectrum 
arises, when there is a symmetry which keeps together 
the Fermi points of different fermionic species, i.e. pro- 
hibits splitting into elementary Fermi points, i.e. Dirac 
points with N — 1. Similar symmetries exist between 
the fermions in the Standard Model of particle physics, 
see [31] and references therein. 

There is a major difference between our approach 
and that of Kitaev who studied the classes of Hamil- 
tonians [TU]. Instead of classification of Hamiltonians 
we make the topological classification of the ground 
states taking into account the symmetry of the ground 
state. We use the top-bottom procedure [33], i-e. con- 
sider from the very beginning the microscopic interact- 
ing many-body system, and consider its properties at 
low energy. The objects of the topological classifica- 
tion are response functions in the background of the 
ground state, which belongs to some symmetry class. 
We use the single particle Green's function, but the 
higher order Green's function can also be appropriate. 
The single particle Green's function matrix at zero en- 
ergy represents the effective single-particle Hamiltonian 
H(p) = G _1 (cj = 0, p), which by definition is non- 
interacting. In Refs. [351 [371 135] the opposite procedure 
is used: the authors start with non-interacting Hamilto- 
nians and extend them to the interacting Hamiltonians. 



This procedure makes sense when one constructs the 
artificial fermionic or bosonic systems. But it is not ap- 
plicable to the natural many-body systems which are 
interacting on the microscopic level: such bottom-top 
procedure is not unique, and it is unable to restore the 
original microscopic Hamiltonian. 

Our approach is also applicable to the vacuum of 
particle physics. The Lorentz symmetry prohibits the 
splitting of the Dirac points, but it also prohibits the 
non-linear non-relativistic spectrum. Situation changes 
if the Lorentz symmetry is viewed as an emergent phe- 
nomenon, which arises near the Dirac point (Fermi 
points with N — ±1). In this case both variants are pos- 
sible, depending on symmetry: splitting of Dirac points 
or formation of non-linear non-relativistic spectrum in 
the vicinity of the multiple Fermi point. In both cases 
the mixing of fermions violates the effective Lorentz 
symmetry in the low-energy corner. This phenomenon, 
called the reentrant violation of special relativity |22j . 
has been discussed for Np = 3 fermion families in rela- 
tion to neutrino oscillations [35] • Influence of possible 
discrete flavor symmetries on neutrino mixing has been 
reviewed in Ref. [40] . 
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